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A fundamental element of quantum information processing with photonic qubits is the nonclassical quantum 
interference between two photons when they bunch together via the Hong-Ou-Mandel (HOM) effect. Ulti¬ 
mately, many such pure photons must be processed in complex interferometric networks, and for this it is 
essential to synchronize the arrival times of the flying photons preserving their purity. Here we demonstrate for 
the first time the HOM interference of two heralded, pure optical photons synchronized through two indepen¬ 
dent quantum memories. Controlled storage times up to 1.8 /xs for about 90 events per second were achieved 
with purities sufficiently high for a negative Wigner function confirmed with homodyne measurements. 


Optical photons are a fundamental resource to encode flying 
quantum bits for quantum communication and computation. 
In particular, in linear-optics quantum information processing 
m, universal two-qubit gates rely upon nonclassical quan¬ 
tum interferences, where photons tend to bunch due to their 
bosonic nature. The elementary manifestation for this is the 
so-called Hong-Ou-Mandel (HOM) effect Q: when two in¬ 
distinguishable single photons |1,1) enter a balanced beam 
splitter, they bunch in either of the two output ports, resulting 
in a HOM state |HOM(6>)) = (|2,0) - e2*<^|0,2))/V2 with 
some relative phase 0. For large-scale quantum computation, 
many pure single photons must be available simultaneously at 
the input ports of large interferometric networks, in order to 
apply the corresponding gate sequences at the same time on all 
the initial qubits. Numerous tests of the HOM effect have been 
performed over the past few decades mainly in order to char¬ 
acterize single-photon sources, such as parametric down con¬ 
verters O, trapped single neutral atoms fH, ions lO, atomic 
ensembles i), quantum dots 13, and nitrogen vacancy cen¬ 
ters in diamond M- However, the simultaneous occurrence of 
two single photons at a beam splitter has depended on a ran¬ 
dom coincidence between two independent statistical sources. 
One possible way toward scalability is to combine statistical 
photon sources with quantum memories, by which a photon 
is stored until the other photon is available too EM- How¬ 
ever, such quantum storage must not be at the expense of the 
single-photon purity. Only in a very recent experiment, optical 
single-photon states sufficiently pure to show a negative dip 
in their Wigner functions m were released in a controlled 
fashion (quasi on demand) from a memory system based on 
optical cavities |[l2l[T3 or atomic ensembles lIT^ . 

Here we report the next significant step beyond this: the 
controlled HOM interference of two nearly pure photons that 
emerge from two independent quantum memories employing 
the cavity-storage method |[l2l[T3. The resulting states still 
have single-photon purities above 0.5, which is sufficient for 
a negative Wigner function. Thanks to the memories, utiliz¬ 
ing controlled storage times of up to 1.8 /iS, the output HOM 
state can be synchronized. Compared to previous works with 


atomic memories, the memory times of our all-optical system 
are of similar order miioiiia, however, the purities of the 
synchronized photons are raised to an unprecedented, qual¬ 
itatively different level. We believe that this controlled, al¬ 
most on-demand demonstration of the HOM effect represents 
a breakthrough toward scaling up photonic quantum inter¬ 
ference experiments, with direct applications in linear-optics 
quantum computation (H, quantum communication ca, and 
boson sampling C3. 

Another important aspect of our approach, making it dis¬ 
tinct from all postselective schemes based on particle-like 
click-by-click photon detections, is the characterization of the 
resulting HOM state from a wave-like perspective using ho¬ 
modyne measurements of field quadrature amplitudes. Here 
we will show that a characteristic pattern in the wave basis 
survives even after the active synchronization with quantum 
memories. A fundamental feature of quantum mechanics is 
the wave-particle duality, and our demonstration looks at the 
famous HOM effect from a completely wave-like angle. This 
is not only of fundamental interest but also practically im¬ 
portant for optical hybrid quantum information processing, 
where both continuous wave and discrete particle properties 
are exploited for quantum state preparation, processing, and 
detection ESI 091. Since the first demonstration El, the HOM 
effect has been always demonstrated as a dip in the coinci¬ 
dence probability of photon detections at both output ports 
of the beam splitter. However, the HOM dip only reflects a 
particle-like aspect of the HOM state, in which photonic par¬ 
ticle bunching becomes manifest. The more general, actual 
quantum nature of that prominent optical quantum state, such 
as quantum entanglement, cannot be revealed by correlation 
measurements in a fixed particle basis. The wave-basis im¬ 
age of the HOM state that we present here as a counterpart 
to the particle-basis HOM dip is a correlation pattern similar 
to a four-leaf clover that vanishes and reappears depending 
on the relative phase 6»Eq1|2D. The phase-dependent pattern 
is actually sufficient for fully characterizing the HOM state 
and obtaining its density matrix (T2 \ . However, the pattern is 
very fragile against optical losses and detection noises mostly 


because the homodyne detection is sensitive to vacuum fluc¬ 
tuations, unlike the HOM dip whose shape is in principle un¬ 
changed even for large optical losses. Therefore, in order to 
observe the phase-sensitive clover pattern, highly pure single 
photons must be prepared simultaneously and detected with 
very low-noise homodyne detectors, which became possible 
only very recently |[T2Hl4l[20l|2Tl. 

Our experimental setup is schematically shown in Fig. la, 
where every one of the two single-photon sources is enclosed 
by a memory system (Memory-1 or Memory-2) which is com¬ 
posed of two concatenated cavities |[l2l[T3. Each individual 
single-photon creation corresponds to an ordinary quantum 
optical heralding scheme, where photons are probabilistically 
but simultaneously produced in pairs by nonlinear optical ef¬ 
fects and one photon serves as the herald of the other ITTIl . 
Single-photon sources based on such nonlinear optical effects 
have good controllability of the wavelength unlike other on- 
demand-type sources. The special interference inside the con¬ 
catenated cavities ensures that the photon to be measured is 
released to the outside while the photon to be prepared stays 
inside. After heralding, the stored single photon is released 
on demand by rapidly switching the cavity resonance via an 
electro-optic effect flJl [TSl, which is also different from a 
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FIG. 1. Experimental setup and characterization, (a) Experimental 
setup. EOM, electro-optic modulator; APD, avalanche photodiode; 
LO, local oscillator; RM, replaceable mirror, (b) Experimentally 
estimated wavepackets of single photons released from Memory-1 
(blue traces) and Memory-2 (red traces) after various fixed storage 
times. The plots are vertically shifted according to the storage times 
from 0 ns to 500 ns by 50-ns steps, respectively. The horizontal axis 
is the time relative to each heralding event. The intrinsic delay in 
the memory release is about 50 ns. (c) Experimental single-photon 
purities after the synchronization without any correction of detec¬ 
tion inefficiencies. Circles and squares are (l|Trfc{B^pout^}|l) for 
Memory-1 (k = 2) and Memory-2 (k = 1), respectively. They are 
calculated from the output with numerical inverse of the balanced 
beam splitter B. 


simple storage-loop switching scheme 1^ . 

In a preliminary experiment, we tested the performance of 
the individual memory systems by using a highly reflective 
mirror as the replaceable mirror (RM) in Fig. la. For both 
input ports, we estimated the wavepackets of the released sin¬ 
gle photons for various flxed storage times after the herald¬ 
ing (TJl [TSl [24]|. Figure lb shows the longitudinal modes 
of the released wavepackets in the time domain, vertically 
shifted depending on the corresponding storage times. It can 
be confirmed that the single-photon wavepackets are correctly 
shifted by the memories without deformation, and also that 
the wavepackets from the two independent memories are al¬ 
most identical. This sameness (indistinguishability) of the 
wavepackets is critical for the HOM interference effect. 

In the HOM experiment, we synchronized the photon 
sources and then obtained an output HOM state. We acquired 
two-mode quadrature data for 390,636 events when both sin¬ 
gle photons were heralded within 2 /is. This large number of 
events was acquired within only 3 hours. The storage time 
of each event was recorded together with the homodyne out¬ 
comes. Then, in order to observe the dependence on the max¬ 
imum storage time r^ax, we analyzed those events associated 
with a storage time between 0 and Tmax- Figure Ic shows 
the individual input single-photon purities, calculated with the 
maximum-likelihood method Ga for each Tmax- Although the 
single-photon purities degrade for longer storage times due to 
a finite memory lifetime ca, the purities were kept above the 
0.5-bound of the Wigner-function-negativity condition for up 
to about 1.8 /is, which corresponds to about 90 events per sec¬ 
ond. The coherence time of the wavepackets for the 0.5 crite- 
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FIG. 2. Theoretical two-mode quadrature distributions in the ideal 
cases, (a) |2,0) state, (b) |0,2) state, (c and d) Coherent superpo¬ 
sition |HOM(^)) for ^ = 0° and for 0 = 90°, respectively. The 
probability amplitudes at the origin interfere destructively (c) or con¬ 
structively (d). The marginal distributions are also shown in the top 
and right panels. 


2 






































































ria without synchronization is estimated at 72 ns ca ,and thus 
the memory enhancement in the event rate can be calculated 
as a factor of 25. 

Let us now discuss the image of the HOM state in the wave 
basis, as shown in Fig. 2. Theoretically, vacuum |0), single¬ 
photon |1), and two-photon |2) states have one, two, and three 
peaks in their quadrature distributions, respectively, and thus 
the two-dimensional distributions of the separable states |2, 0) 
and |0, 2) look like those in Fig. 2a and Fig. 2b, respectively. 
When 12,0) and |0,2) are coherently superimposed to ob¬ 
tain the HOM state, |HOM((9)) = (|2,0) - 2 ))/a/2, 

the probability amplitudes at the origin destructively or con¬ 
structively interfere depending on 0 Il20l|2l]|. As a result, 
the remaining distribution patterns are the four-leaf clover 
for 6> = 0° (Fig. 2c) and a concentric pattern for 0 = 90° 
(Fig. 2d). 

These characteristic distribution patterns are indeed ob¬ 
served in our experiment by homodyne detections. In 
Figs. 3a,b, the distributions of the output for 0 = 0° and 
90° are shown, together with the individual histograms of the 
single-mode quadratures. These are the results for a r^ax of 
400 ns as a typical example. The output barely exhibits side 


fringes reflecting a large two-photon component (top and right 
panels in Figs. 3a,b). The two-mode distributions of the out¬ 
put completely change depending on 0. The four-leaf clover is 
most pronounced at = 0° whereas the clover totally disap¬ 
pears at 6> = 90° as expected from Figs. 2c,d ||20l|2Tl. These 
distributions are totally different from that of the case without 
synchronization (Fig. 3c). The results here are only particu¬ 
lar examples, and the full results are presented in the supple¬ 
mentary information CD, together with those for the phase- 
independent input dual-single-photon state. 

So far, we have discussed the relative phase 0 of the HOM 
state |HOM(6»)) = C/(6»)|HOM(0)), where U{e) = jg 

the phase-shift operator acting on the second mode. From a 
different viewpoint, this can be reinterpreted as a phase shift 
of the measured quadratures {xi^X 2 \U{ 0 ) for the fixed state 
|HOM(0)), since Pr(a;i, 2 : 2 ) = |(a;i,a;2|C/(6»)|HOM(0))|l In 
the latter interpretation, the wave-basis two-mode distribu¬ 
tions with various measurement phases contain the complete 
information for estimating the quantum state 1^ . In Fig. 3d, 
the density matrix of the output state pout is shown in the num¬ 
ber basis. The ideal density operator of the HOM state is 
Phom = |HOM(0))(HOM(0)| = i(|0,2)(0,2| + |2,0)(2,0|- 
|0, 2) (2, 0| — |2, 0) (0, 2|), and this structure also becomes 
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FIG. 3. Experimental outputs, (a and b) Experimental quadrature distributions of the output state after the synchronization for ^ = 0° and for 
0 — 90°. (C) Experimental distribution without synchronization for ^ = 0°. (d) Density matrix of the output state after synchronization. The 
imaginary part of the density matrix is omitted because it is negligibly small; it can be found in the supplementary materials Ca. (e) Cross 
section of the output Wigner function, (f) Density matrix of the output without synchronization. 
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apparent in the experimental output state pout- the diago¬ 
nal elements of |2,0)(2,0| and |0,2)(0,2| indicate the pho¬ 
ton bunching effect on the |1,1)(1,1|-component of the in¬ 
put ca, while the presence of the off-diagonal, negative el¬ 
ements of |0, 2) (2, 0| and |2, 0)(0, 2| proves that the bunched 
components are quantum-mechanically superimposed rather 
than classically, incoherently mixed. These off-diagonal el¬ 
ements are (necessary and) sufficient for the entanglement 
between the two output modes. Quantitatively |[25l . a loga¬ 
rithmic negativity log 2 ||pout^^|| = 0.37 ± 0.01 is obtained, 
while log 2 IIphom^HI = 1 in the ideal case. Moreover, 
strong nonclassicality of the output state is evident from its 
negative Wigner function lT^(xi, pi, X 2 , P 2 ) • The cross sec¬ 
tion iy(xi, 0, X 2 , 0) shown in Fig. 3e has negative values in 
a structure reflecting the four-leaf-clover pattern. A further 
evaluation can be found in the supplementary information 

ca. 

In conclusion, we have experimentally demonstrated the 
famous HOM interference effect in a way that is potentially 
scalable for large-scale quantum information processing. To 
achieve this, two nearly pure single photons enter a beam 
splitter almost on demand after their synchronized release 
from two independent optical (cavity-based) quantum mem¬ 
ories. We have also shown the wave-like aspect of the re¬ 
sulting HOM state via homodyne detections, characterized 
by a phase-dependent appearance of a four-leaf clover pat¬ 
tern. The quality of the states is sufficient to exhibit a nega¬ 
tive Wigner function and entanglement even after the memory 
storage. The HOM state interpreted as a NOON state 1261 is 
also potentially useful in quantum metrology and sensing. 
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SUPPLEMENTARY METHODS 
Details of experimental setup 

The light source of our experiment is a continuous-wave 
(CW) Ti:sapphire laser operating at the wavelength of 860 nm. 
In addition to the setup shown in Fig. la, there are two opti¬ 
cal cavities which are omitted from Fig. la. One cavity is a 
second-harmonic generator, which is a bow-tie-shaped cavity 
and contains a periodically-poled KTi 0 P 04 (PPKTP) crystal 
as a nonlinear optical medium. The resulting continuous out¬ 
put beam at the wavelength of 430 nm is, after a frequency 
shift by an acousto-optic modulator, directed as a pump beam 
to each memory cavity, which contains a periodically-poled 
KTi 0 P 04 (PPKTP) crystal as a nonlinear optical medium 
and works as a non-degenerate optical parametric oscilla¬ 
tor (NOPO). The pumping power at each memory cavity is 
about 3 mW. Optionally, the pump beams can be individually 
blocked after a herald in order to prevent further production 
of photons, but this was not employed in this demonstration, 
because the probability of such unwanted events is small. The 
other cavity is a mode-cleaning cavity, by which the transverse 
mode of the local oscillators is purified to a TEMqo mode in 
order to maximize the interference visibility at the homodyne 
detections. The power of each local oscillator is about 18 mW. 

The concatenated-cavity-based memory systems (Memory- 
1, Memory-2) have the same design, similar to our previous 
experiment ca but with parameters slightly different from 
before. The design of the concatenated-cavity-based memory 
system is as follows. The transmissivity of the coupling mir¬ 
ror that couples the memory cavity and the shutter cavity is 
1.6 % and the transmissivity of the outcoupling mirror at the 
exit of the shutter cavity is 24.7%. The former transmissiv¬ 
ity corresponds to a compromise between long lifetime and 
phase-locking stability, while the latter is set in accordance 
with the former in order to operate the system nearly at the 
critical damping condition Cll. The memory cavity has a free 
spectrum range (FSR) of 214.1 MHz. Spontaneous paramet¬ 
ric down conversion inside the memory cavity produces signal 
and idler photons which are separated by this FSR. The shut¬ 
ter cavity contains an electro-optic modulator (EOM), which 
is a RbTi 0 P 04 (RTF) crystal with an aperture size of 4 mm x 
4 mm. It is driven by a high-voltage switch (Bergmann Mess- 
gerate Entwicklung KG), whose voltage is around 900 V to 
match the frequency shift to the FSR of the memory cavity. 


In order to stabilize the whole setup, all of the resonant 
frequencies of the cavities and the phases of the local oscil¬ 
lators are each electronically controlled by using a feedback 
loop. Each controller is composed of an analog feedback 
controller, an error detector of locking, and a digital scanner 
for the error recovery. The error detection signals from all 
of the controllers are brought together by logical OR gates, 
and then sent to the timing controller in order to pause the 
measurement when there is an error. The timing controller is 
a field-programmable gate array (FPGA) (Virtex-4, Xilinx), 
which processes the heralding signals and controls the timing 
of the photon release. The same timing controller actually also 
controls the sequence of the switching between the feedback 
phase and the measurement phase, as explained below. For the 
analog feedback control, we monitor the optical systems by 
means of bright beams. However, such bright beams, except 
for the local oscillator beams and the pump beams, represent 
an extra complication in our single-photon-level experiment. 
Therefore we cyclically switch the optical systems. One phase 
is the feedback phase where the bright beams are injected to 
the cavity systems. The other is the measurement phase where 
the bright beams are blocked and the two-photon interference 
is tested. The switching rate is 5 kHz, and the duty cycle is 
40% for the measurement phase. 

Timing control sequence 

As mentioned above, the FPGA controls the release tim¬ 
ing of photons. The FPGA clock frequency is 100 MHz, and 
thus, in our system, the release timings are synchronized with 
10 ns intervals. Note that the timing jitter of the driving sig¬ 
nals in the timing controller was negligibly small (less than 
1 ns), compared to the width of the wavepackets of photons 
(about 100 ns). The photon arrival times of the two inputs are 
matched by adjusting the electric cable lengths from the tim¬ 
ing controller to the EOMs, as well as the optical path lengths 
from the exits of the shutter cavities to the balanced beam 
splitter of the HOM interferometer. The lengths of the two 
output arms to the data-storage oscilloscope are also matched. 
These were tested in the preliminary experiment. 

In the HOM experiment, the photons are released when 
both heralding events of Memory-1 and Memory-2 happen 
within 2 /i5. If the second herald does not occur by 2 /is, 
the timing controller becomes idle and waits for 5 /iS, which 
is long enough compared to the memory lifetime of about 2 


5 


/iS. During this dead time, the heralded photon is almost lost 
and the memory systems are almost initialized again. 


SUPPLEMENTARY DISCUSSION 


Detailed analysis of the preliminary experiment 


In the preliminary experiment, we tested our memory sys¬ 
tems by using a highly reflective mirror as the replaceable 
mirror (RM). For each Memory-/c, we repeated single-photon 
generation and homodyne detection 43,404 times for each 
variable set storage time r after the herald, changed from 0 
ns to 500 ns at 50-ns intervals. From the results, we estimated 
the shapes of the wavepackets by utilizing principal 

component analysis Ha [m Ian, as shown in Fig. lb, and cal¬ 
culated the single-photon purities P/c(t) (corresponding to the 
respective single-photon fractions), as shown in Fig. SI a. The 
single-photon purities are shown to degrade for longer stor¬ 
age times indicating a flnite memory lifetime, but the Wigner- 
function-negativity condition ifTTlfT^ with purities above 0.5 
is still satisfled for up to about 0.4 /is. As shown in Fig. Ic, the 
maximum storage time r^ax of the synchronization can be set 
much longer than this 0.4 /is, thus preserving the negativity. 

The shapes of the estimated wavepackets //c(t;r), shown 
in Fig. lb, are independent of r, excluding pre-leakage. Fur¬ 
thermore, they are almost the same between Memory-1 and 
Memory-2. The overlap C(0) = | J 0 )/ 2 (t; 0)p was 

0.992. We deflned the signal wavepacket for the HOM exper¬ 
iment as an average temporal mode: fuouit) = 0) + 

/ 2 (f; f 0) + / 2 (t; 0)p. The overlap between 


the signal wavepacket and the estimated wavepacket released 
from the individual Memory-/;: after r, deflned as C/c(r) = 

I f dtfuoM (f ~ //c (f; ^) P 5 was larger than 0.96 for up to r = 
500 ns, as listed in Supplementary Table ^ This high value 
contributes to the visibility of the HOM interference and the 
purities of the input states. 


The estimated purities P/c(t) corresponding to the esti¬ 


mated wavepackets fk{t; r) in Fig. lb are listed in Table SI 


corresponding to the exact values of the points in Fig. SI a. In 
order to exclude the influence of the mode mismatches on the 
input purities, we also calculated purities by using the 

time-shifted signal wavepacket fuouit — r), as also shown in 
Table From the values below, we can see that the con¬ 
tribution of the pre-leakage is very small. The decay of the 
purity is thus due to the optical losses in each memory sys¬ 
tem. These results are fitted with an exponential decay curve 
Pk{^) = A:(0) exp(—r/f/c), where P/c(0) means the initial 
purity and means the memory lifetime for the memory sys¬ 
tem. 


Table SI. Overlap of the wavepackets and purity for each 
wavepacket. 


Set storage time 

Overlap 

Purity (fk{t;T)) 

Purity (fnouit - r)) 

T (ns) 

Ci{r) 

C2{t) 

Pi{r) 

P2{t) 

Plir) 

Piir) 

0 

0.998 

0.998 

0.606 

0.642 

0.608 

0.642 

50 

0.994 

0.992 

0.584 

0.613 

0.584 

0.613 

100 

0.994 

0.987 

0.579 

0.596 

0.577 

0.593 

150 
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Fig. SI. Input characterization and dependence on the 
maximum storage time, (a) Input single-photon purities 
(l|p|l) with respect to the storage time. Blue open circles 
and red open squares represent Memory-1 and Memory-2. 

Error bars are ±0.006. Traces (i and ii) are exponential 
fittings, (b) Number of events for the input state (blue) and 
the output HOM state (red) with respect to the maximum 
storage time Tmax- 


In the HOM experiment, we synchronized the photon 
sources, created an initial dual-single-photon state, and then 
obtained an output HOM state by replacing the highly- 
reflective RM in Fig. la with a balanced beam splitter. The 
storage time of each event was recorded together with the ho¬ 
modyne outcomes. Then, in order to observe the dependence 
on Tniax, we analyzed those events associated with a storage 
time between 0 and Tmax- Figure Sib shows the obtained num¬ 
ber of events with respect to each Tmax- It indicates a linear 
increase of the dual-heralding event rate for any extension of 

±nax* 

Here we discuss how the maximum storage time Tmax is de¬ 
pendent on the dual-heralding event rate. We assume that the 
individual heralding events of the Memories occur at random 
and independently, i.e. they follow the Poisson distributions 
with average event rates Rk for each Memory-Zc (/c = 1, 2). 
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Considering the (dimensionless) measurement duty S and ig¬ 
noring the 5-jlis dead time for simplicity, the theoretical dual- 
heralding event rate is given by i?duai(^max) = 
when Tmax is much smaller than S/Ri and Here the 

factor of 2 effectively stems from the fact that the first herald¬ 
ing event may occur in either of Memory-1 and Memory-2. 
This linear relationship of the dual-heralding event rate for 
^max can be confirmed in Fig. Sib. The experimental pa¬ 
rameters were i ?2 = 3, 200 counts per second (cps), and 
S = 40%. From above, the predicted dual-heralding event rate 
for Tmax = 1.8 /i 5 is i?duai(l -8 = 92 cps, which agrees 
well with the experimental value of about 90 cps. 

Coherence time of wavepackets 

In order to estimate the enhancement by the memories in 
the dual-heralding event rate, here we discuss the coherence 
time of the wavepacket without synchronization, from the as¬ 
pect of the 0.5-bound of the Wigner-function-negativity con¬ 
dition. However, we note that this is a problem dependent on 
how the signal wavepacket is defined. A reasonable definition 
of the signal wavepacket is to set the temporal origin of the 
signal wavepaket to the mean of the two heralding event tim¬ 
ings: when individual heralding events happen at ti = 0 and 
t 2 = At, then a signal wavepacket /hom(^ — At/ 2 ) has a bal¬ 
anced overlap with both of fi (t; 0) and /2 (t — At; 0). By us¬ 
ing this definition, we calculate the average purity for various 
timing-mismatched events (—Tcoh < At < Tcoh) as a function 
of the maximum time difference Tcoh- Then we find that the 
average purity crosses the 0.5-bound of the Wigner-function- 
negativity condition at Tcoh = 22 ns. Comparing this Tcoh of 
72 ns with the Tmax of 1.8 /is for the synchronized case, and 
taking into account the linear dependence on r^ax, the mem¬ 
ory enhancement in the event rate can be estimated at a factor 
of 25. 


Two-mode wave functions and probability distributions 

Let us now discuss the image of the HOM state in the 
wave basis by considering wavefunctions (probability ampli¬ 
tudes) in Fig. S2. Theoretically, vacuum |0), single-photon 
| 1 ), and two-photon | 2 ) states have zero, one, and two nodes 
in their Schrodinger standing wave functions 'ipn{x) = {x\n) 
for the field quadrature amplitude x, respectively (Fig. S2a). 
The two-mode wave functions of separable states |2,0), |0,2), 
and | 1 , 1 ) are products of the single-mode wave functions 
(Figs. S2b-d). When |2,0) and |0,2) are coherently super¬ 
imposed to obtain the HOM state, |HOM(^)) = (|2,0) — 
2 ))/\/ 2 , the probability amplitudes at the origin de¬ 
structively or constructively interfere depending on 0. As a 
result, the superimposed wave function corresponds to a four- 
leaf clover for 0 = 0° (Fig. S2e) and a concentric pattern for 
0 = 90° (Fig. S2f). Comparing Fig. S2e with Fig. S2d repre¬ 
senting the input state 11 , 1 ), we can see that the interference at 



Fig. S2. Wave functions and probability distributions, (a) 
Single-mode wave functions V^n(^) of photon-number states, 
(b to f) Two-mode wave functions. Those of separable states 
| 2 , 0 ), | 0 , 2 ), and input | 1 , 1 ) are products of single-mode 
wavefunctions {xi, X 2 \ni,n 2 ) = (^i)' 0 n 2 (^ 2 )- At the 

HOM output, the probability amplitudes at the origin 
interfere destructively (e) or constructively (f). (g to i) 
Probability distributions Pr(xi, X 2 ), which are absolute 
squares of the wave functions (d to f), respectively. The 
marginal distributions are also shown in the top and right 
panels. 


the balanced beam splitter for 0 = 0° corresponds to a 45° ro¬ 
tation of the pattern. The probability distributions Pr(xi, X 2 ) 
are the absolute squares of the wave functions (Figs. S2g-h). 


Experimental quadrature distributions for the 36 measurement 

bases 

In the HOM experiment, the optical phases ( 6 >i, 6 ^ 2 ) of the 
two measured quadratures {xi ( 6 >i), X 2 { 02 )) are independently 
set to 0°, 30°, 60°, 90°, 120°, and 150°. The quadrature dis¬ 
tributions for the tested 6 x 6 = 36 combinations are shown in 
Fig. S3. These are the results for a Tmax of 400 ns as a typical 
example. The distributions for equal relative phases, as can 
be seen in the plots along the diagonal of Fig. S3, are almost 
the same, thus refiecting the phase-insensitivity of the input 
single-photon states. The characteristic results of the output 
for the relative phase 0 = 0i — O 2 of 0° and 90° are collected 
in Fig. S4, together with those of phase-insensitive input state. 
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Fig. S3. Experimental two-mode quadrature distributions of the output state. Quadrature distributions for 36 combinations of 

the optical phases (Oi, O 2 ) are shown. 
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Fig. S4. Experimental quadrature distributions of the input and output states, (a to c) Experimental quadrature distributions of 

the input state and the output state for 6> = 0° and for 0 = 90°. 


Quantitative evaluation of the entanglement of the output HOM 

state 


entanglement. The intensity cross-correlation factor, 
which corresponds to the depth of the HOM dip, 
(a|a 2 a 2 ai)/(a|ai)(a 2 a 2 ), is 0.24 ± 0.01. The interference 
visibility, defined as (|a|ai — a 2 a 2 |)/(a|ai + a 2 a 2 ) is 0.902 
± 0.005. This is much higher than the phase-randomized co- 


Here we quantitatively evaluate the output state with 
regards to its correlations, interference visibility, and 
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herent light limit of 0.5, which conhrms the non-classicality 
of the output state. In order to verify the quantum entangle¬ 
ment in the two output ports, we calculated the entanglement 
criterion Ell. The logarithmic negativity 1^ of the output 
state was log 2 | |pout^^ || =0.37 ±0.01, where the superscript 
of Ti means partial transpose with respect to mode 1. In or¬ 
der to exclude the contribution of undesired photon-number 
components to the entanglement, we also calculated the log¬ 
arithmic negativity of the output state after numerically ap¬ 
plying local hlters |0)(0| ± |2)(2| to individual modes. The 
local hlters do not increase the entanglement but just remove 
irrelevant photon-number components and transform the out¬ 
put state to = PsPs ± (1 — Here, ps is the hl- 

tered density operator of pout in the Hilbert space spanned by 
|0, 0), |0, 2), |2, 0), and |2, 2), is the probability fraction of 
the hltered components which was experimentally 0.470, and 
p^ is a non-entangled state in the orthogonal Hilbert space. 


For this hltered output state, the logarithmic negativity was 
iog 2 IIPout^MI = 0.30 ± 0.01. From comparison of these 
two values, we can see the main contribution of the entan¬ 
glement came from the superposition of |2, 0) and |0, 2), but 
there is also a small discrepancy of the log-negativities which 
is mainly due to the higher-order photon components |2,1) 
and |1, 2) of the input states. The non-zero log-negativities of 
the hltered and unhltered states conhrm the entanglement of 
our output HOM state. 

Simulation of our HOM experiment from the input states 

By using the density matrices of the input states 
(Fig. S5 a,b), we theoretically calculated a density matrix 
of the output state with the ideal balanced beam splitter 
(Fig. S5 e,f). Compared to this simulation, the density matrix 
of the experimental output HOM state has a similar structure 
(Fig. S5 c,d). 



Fig. S5. Real and imaginary parts of the density matrices, (a, b) Experimental input state, (c, d) Experimental output state, 
(e, f) Simulated result of the output state from the experimental input state with a numerical balanced beam splitter. 
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Wigner function of the output state 


Additional important information about the phase-space 
characteristics of the entangled HOM state can be obtained 


by analyzing its Wigner function. We calculated the Wigner 
function of the output state in order to verify the strong non- 
classicality of this state. The Wigner function of a two-mode 
state p is defined (h = 1) in four-dimensional space as 


W{xi,pi,X2,P2) = jj dyidy2e^ 


{Piyi+P2y2) 


yi y2 


, yi , 2/2 

Xl + -,X2 + - 


where Xk = Xk{0k = 0),pk = Xk{0k = 90°), {k = 1,2) 
correspond to a pair of rotated quadratures. In Fig. S6, the 
Wigner function of the ideal HOM state and that of the exper¬ 
imental output state are shown, for different vertical scales. 
As the Wigner function is very sensitive to the non-classical 
properties of a quantum state, its value quickly changes ow¬ 
ing to experimental imperfections. However, the Wigner func¬ 


tion of the output state still shows the distinctive shape of the 
ideal HOM state: it has a positive value at the origin and 
negative values at four points around the origin in Fig. S6c, 
which come from the strongly non-classical single-photon in¬ 
puts, thus confirming that the output state maintains the strong 
non-classicality after the photon interference. 






Fig. S6. Cross sections of the two-mode Wigner functions, (a, b) Ideal HOM state Phom for comparison, (c, d) The 
experimentally reconstructed output state for Tmax = 400 ns. Note that the vertical scales in (a, b) are different in (c, d). The 
fixed cutting plane {jpi = 0,p2 = 0) is chosen to make the quadrature correlations visible (a, c). The fixed cutting plane 

{jpi = 0, X 2 = 0) is chosen (b, d). 
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